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0.  Introduction 

Mi  . . .  one  half  of  the  Spring  1983  semester  at  ILe  Institute 

* 

for  Mathematics  and  its  Applications.  During  that  time  he  interacted  with 
colleagues,  engaged  in  research,  and  gave  two  public  lectures  at  the  Institute: 
9  (ty  Chaos  in  Phase  Transitions,  a*.  ^ 

9m-  Dynamics  pf  Phase  Transitions. 

- —  ^  ^ 

The  main  thrust  of'Slcmrod^s  research  was  in  two  areas,  specifically: 


Deterministic  chaos  in  materials  exhibiting  phase  transitions,  and 

"  *  r 

$2-411).  Admissibility  criteria  for  weak  solutions  of  the  non -hyperbolic 
conservation  laws  which  describe  dynamic  phase  transitions. 

A  detailed  description  of  this  research  is  given  below.  \ 


AIR  FORCE  0F7I 
NOTICE  OF  TK.V.V 
This  teMirlc  ' : 
approved  for ' 
Distrib\:t,  1 

VATTOEW  J.  KEK 
Chief,  Teel r<1: 


C £  i>F  SCIENTIFIC  RESEARCH  l  A* SO  ) 

C”  T  TTAL  TO  DTIC 

’•  -p'> \~t.  l;a«  t'oci  reviewed  tad  IS 

:c  If -V  A >'R  190— 12. 

■  .  1  n:  •- 1  ■'  .1 . 

e  '.l  1  -iforaiat  1  on  division 


i 


1 


1.  Chaos  In  Materials  Exhibiting  Phase  Transitions 

Id  recent  work  Holmes  and  Marsden  (1981)  have  employed  the  Mel'nikov 
(1963)  technique  to  establish  the  existence  of  deterministic  chaos  in  periodically 
forced  evolution  equations.  In  this  work  Slemrod  (in  collaboration  with  J.E. 
Marsden  of  the  University  of  California,  Berkeley)  combined  this  program  by 
analyzing  two  problems  arising  from  the  vander  Waals  theory  of  phase  transi¬ 
tions  (see  vender  Waals  (1893,  1979)]. 

The  first  problem  considered  was  dynamic  spinoidal  decomposition. 
Assume  we  have  a  vander  Waals  fluid  with  isotherms  as  shown  in  Figure  1. 


Here  p  denotes  the  pressure  given  by  the  vander  Waals  constitutive 
relation 


p(w,9) 


R9  _a_ 

w-b  "  2  • 

w 


(1.1) 


R,  b,  a  are  positive  constants,  p  is  the  density,  w  ■  1/p  the  specific  volume, 
and  9  is  the  absolute  temperature.  Figure  1  sketches  the  vander  Waals  isotherms 
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for  0  above,  equal  to,  and  below  the  critical  temperature  3crit  *  8a/27bR. 

[A  good  reference  Is  Fermi  (1936J.3  For  0^  <  0  ,  p(w,0)  has  the  following 

features: 

(0  pw(w,eQ)  <0  on  (b,  o)U(/3.  C d)  , 

(ii)  pw(w,0Q)  >0  on  (a,  0)  , 

(iii)  Pw(a,  0Q)  -  PwO,  0Q)  *  0  . 

The  domain  (b,o)  corresponds  to  the  fluid  being  liquid;  the  domain  O.oo) 

corresponds  to  the  fluid  being  vapor;  the  domain  (a,  Pi  is  the  unstable  region  and 

is  referred  to  a6  the  spinoidal.  Also  the  point  w  is  where  p  =  0  on  the 

U  4  3^ 

graph  of  the  0^  isotherm  (w^  is  the  zero  of  R0qW  +  a(w  -b)  =  0) . 

The  fluid  flow  is  thought  of  as  taking  place  along  the  x-axis  in  a  tube  of 

unit  cross  section  of  fixed  volume.  If  0,  >  0  the  state  w  =  w.  will  describe 

1  cnt  0 

a  stable  homogeneous  configuration  for  the  hydrodynamic  equations  [see  Felder- 
hof  (1970),  Slemrod  (1982)3 .  We  then  "instaneously"  quench  the  fluid  by  reducing 
the  temperature  to  0Q  <  0  .  The  homogeneous  state  wQ  will  now  be  in  the 

unstable  6pinoidal  region  of  the  0^  isotherm.  We  studied  the  effect  of  a  small 
periodic  in  time  fluctuation  of  the  absolute  temperature  about  0^.  Specifically 
we  showed  how  the  loss  of  stability  of  Wq  is  accompanied  by  deterministic  chaos 
in  that  there  will  be  a  Poincare-Birkhoff-Smale  horseshoe  in  the  dynamics. 

The  second  problem  considered  is  the  effect  of  a  small  thermal  perturbation, 
periodic  In  space,  of  the  form  0(x)  *  0^+  €  cosqx,  e  small,  on  the  equilibrium 
configuration  of  an  infinite  tube  of  liquid  under  given  applied  load.  In  thi6  case  we 
show  there  are  solutions  with  the  features  of  both  metastable  and  co-existing 
phases  that  exhibit  spatial  chaos. 

It  seems  interesting  to  note  that  the  above  equilibrium  result  has  features 
in  common  with  experimental  observations  presented  by  Dr.  L.  Zapas  of  NBS  at 
the  IMA  Workshop  on  Orienting  Polymers,  March  21-25,  1983. 

The  results  of  this  research  will  appear  as  an  IMA  report  and  will  be 
submitted  for  publication. 
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2. 

Id  recent  years  an  ever-increasing  number  of  people  have  taken  up  Gelfand' s 
(1963)  program  of  parabolic  regularization  of  quasilinear  conservation  laws.  While 
this  work  has  usually  been  motivated  by  the  search  for  reasonable  admissibility 
criteria  for  weak  solutions  of  the  conservation  laws,  it  has  also  provided  a  tech¬ 
nique  for  proving  existence  of  solutions  as  well.  Specifically  in  this  regard  we 
note  the  work  of  Oleinik  (1957)  for  a  scalar  conservation  law  and  the  recent 
remarkable  results  of  DiPerna  (1983)  for  a  two-dimensional  system  of  conserva¬ 
tion  laws.  In  light  of  this  work  Slemrod  attempted  to  give  (where  possible)  a 
physical  motivation  to  the  Gelfand  program  with  respect  to  the  equations  of  com¬ 
pressible  fluid  flow. 

Consider  a  system  of  n  conservation  laws  written  in  vector  form 

F(U)t  +  G(U)x  -  0  (2. 1) 

where  U(x,  t)  is  an  n-vector,  -co  <  x  <  CD,  and  t  >  0.  As  is  well  known,  the 
initial  value  problem  (2.1)  and 

U(x.O)  -  UQ(x)  (2.2) 

does  not  in  general  possess  global  smooth  solutions  if  (2. 1)  is  nonlinear,  no 
matter  how  smooth  Uq.  Hence  we  are  led  to  search  for  weak  solutions  of  (2, 1). 
Unfortunately,  there  is  a  serious  price  to  pay  for  enlarging  the  solution  class, 
namely  loss  of  uniqueness  of  solutions. 

The  Gelfand  program  of  parabolic  regularization  attempts  to  give  a  way  to 
pick  out  the  physically  relevant  weak  solutions  of  (2. 1)  and  thereby,  hopefully, 
recover  uniqueness  of  weak  solutions  of  (2. 1).  Specifically  Gelfand  suggested  that 
the  physically  admissible  solutions  of  (2.1)  should  be  limits  of  solution  U€  of 
the  "artlfically  viscous"  equations 

F(U€)t+G(U€)x  -  e«Ux)x  (2.3) 

where  D  :R°  — »Fn,  c  >  0  a  small  parameter.  For  example,  the  equations  of 
viscous  isothermal  compressible  fluid  flow  in  Lagrange  coordinates  are 


Admissibility  Criteria  for  Weak  Solutions  of  Conservation  Laws 


(2.4) 


<  +  *e\ 
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Unfortunately  D  in  this  example  is  singular,  whereas  the  recent  results  of 
DiPerna  on  existence  of  weak  solutions  of  (2.1)  as  limits  of  (2.3)  require  D 
nonsingular. 


Slemrod  in  his  research  attempted  to  give  a  physical  justification  for 

nonsingular  D  matrices  based  on  the  vender  Waals-Korteweg  gradient  theory  of 

stress.  U6ing  this  idea  Slemrod  was  able  to  apply  DiPerna’s  technique  to  prove 

the  existence  of  weak  solutions  of  (2.4,  c  *  0)  for  all  t  0  when  p  is  described 

by  the  $  ..  Isotherm  of  Fig.  1. 

J  crit 

This  work  will  appear  as  an  IMA  report  and  will  be  submitted  for 
publication. 
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